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Here we consider the nonlinear evolution of Alfven waves that have been excited by gravitational 
waves from merging binary pulsars. We derive a wave equation for strongly nonlinear and dispersive 
Alfven waves. Due to the weak dispersion of the Alfven waves, significant wave steepening can occur, 
which in turn implies strong harmonic generation. We find that the harmonic generation is saturated 
due to dispersive effects, and use this to estimate the resulting spectrum. Finally we discuss the 
possibility of observing the above process. 
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I. INTRODUCTION 



The launch of large projects for the detection of gravitational waves such as LIGO (Laser Interferometer Gravita- 
tional Wave Observatory) together with ambitious projects under development such as LISA (Laser Interferometer 
Space Antenna) [l| has increased the hope for successful detection of gravitational waves (GW:s) during the next few 
decades, and stimulated much work (e.g. Refs. 0,0)- Naturally the research devoted to detection concerns low am- 
plitude GW:s. In an astrophysical context close to the GW source, the gravitational waves can propagate in a plasma 
medium, and the amplitudes are larger. In Refs. 0, 0, M, B 0> B B E3] the authors have studied nonlinear responses to 
the gravitational wave by the plasma medium, although the backreaction has been neglected. The nonlinear response 
gives raise to effects such as parametric instabilities [!, H, 0, @] j large density fluctuations [!, Q , photon acceleration 
Q and wave c ollap se [Toll . The application of gravitational wave processes to astrophysics has been discussed by for 
example Refs. jl ll. [l2l Tl3L IT3. [Til] . and to cosmology by Refs. [TgI. Il7l. fl8l E^. A number of works studying nonlinear 
propagation of gravitational waves including the backreaction from the plasma have also been written, see e.g. Refs. 

mm _ 

In Ref. [lfj| the coupled evolution of GW:s and (compressional) Alfven waves were considered, and a condition 
for nonlinear wave collapse was derived. Here we develop that work, taking into account the dispersive properties of 
the Alfven waves. For conditions, when wave collapse do not occur, we here show that the nonlinear evolution lead 
to significant wave steepening, and associated high harmonic generation. Furthermore, we find that the harmonic 
generation is saturated due to the dispersive effects associated with the Hall current. As a result, we are able to 
estimate the resulting wave spectrum. Due to the strong harmonic generation, it turns out that electromagnetic wave 
frequencies several orders of magnitudes larger than the original GW frequency can be generated. Finally we discuss 
the possibility that such radiation might be observable with satellite based radio arrays [2l|, [HJ . 



II. COUPLED ALFVEN AND GRAVITATIONAL WAVES 



The metric of a linearized gravitational wave propagating in the z-direction can be written as |23j 

ds 2 = -dt 2 + [1 + h(z - ct)]dx 2 + [1 - h(z - ct)]dy 2 + dz 2 , (1) 

where we have assumed linear polarization, as will be justified below. For an observer comoving with the time 
coordinate, the natural frame for measurements is given by 

eo = d t , e 1 = (l - \h) d x , e 2 = (l + \h) d y , e 3 = d z . (2) 

We now introduce the 3-vector notation such that V = {ei,e.2,e-z), E = (E 1 , E 2 , E 3 ) is the electric field and B = 
{B 1 , B 2 , B 3 ) is the magnetic field. It can be shown [ll[ that in the frame ^ Maxwell's equations can be written 

V E = £- , (3) 
V-B = , (4) 
d t E-c?V-KB = -h. - L , (5) 

dtB + VxE = -^2. . (6) 

ce 
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Here p = J2 S 97" ^ s the charge density, q and n denotes the charge and the proper particle density for a particle of 
species s and j and j are the effective gravitational current densities, defined by 



ceo 



£ 2 



(7) 
(8) 



Next we assume the presence of a background magnetic field, Bo = -Bo e ij and introduce perturbations such that 
n = riQ + 5n, B = (Bq + B x )ei, E = E v e% and v = v y e2 + 1*363, corresponding to compressional Alfven waves 
excited by the GW:s. We also note from Ref. [3( that in the case of gravitational waves propagating in a magnetized 
plasma, with the magnetic field perpendicular to the direction of propagation, only the linear component of the GW 
polarization considered in Eq. {1} couples effectively to the electromagnetic wave. Furthermore v y <C c, (as v y /c ~ h), 
and we therefore neglect terms of the type v 2 and v y h, but allow for v z ~ c. We will also consider slow variations 
such that dt "C w c = qB^/m for each plasma species. 

Following Ref. [Io| . letting the variables depend on z and t, we find that the Maxwell and fluid equations can be 
reduced to 



d t E y - c 2 d z B x + \^ — 7(n + Sn)v y 

d t B x — d z E y 
dt (l(n + Sri)) 
dtilVy) + v z d z (jv y ) 

d t {jv z ) +v z d z (^v z ) 



. '-,1 ,<■)_/, 



~E y d z h+^(B + B x )d t h 
-d z (j(n + 6n)v z ) , 
^-(E v + v z (B Q + B x )) , 

TO 

--» tf (fl0 + B x ) , 



(9) 

(10) 

(11) 
(12) 

(13) 



where Eqs. (|llll - (|13p holds for each particle species. However, for notational convenience we omit the index denoting 
species, as we at this stage want to cover both the case of an electron-ion plasma as well as an electron-positron 
plasma. The above system of Eqs. (P|)- (TT^1) should be complemented by the non-trivial part of the linearized Einstein 
field equations [HUll 



(c~ 2 d 2 t -dl)h = K(5T 11 -5T 22 ) 



(14) 



where n = 8irG/c 4 , G is Newton's gravitational constant and STij are the perturbed components of the energy 
momentum tensor. Provided the coupled Alfven waves and GW:s move together close to the speed of light, the vacuum 
expressions for the effective currents ([7 )) -l[5 |l holds ap pro ximately. The above system of equations has previously been 
studied by Ref. [|| and Ref. [l(| • In particular Ref. jlOl ] derived a system of two coupled equations of the form 



(d t +U(B x )d z )B x 



1 



B d t h , 



-2k- 



BqB x 



/to 



(15) 
(16) 



where U(B X ) = c- (c 3 /2C%) (B / (B Q + 2B X )) 3/2 , C A = c(£ s ^^/^-i/a is the n on-relativistic Alfven velocity and 
uj p = (noq 2 /eoTO,) 1 / 2 is the plasma frequency for each species. In contrast to Refs. flo| we will here concentrate on 
the evolution of the Alfven waves at a distance further from the gravitational source, where the Alfven-GW coupling 
is of less importance, and the GW source terms for the Alfven waves can be n eglected. T hus from now on we will 
assume the perturbations to be of the form ip « ip(z — V^t), where Va — Ca/ v/T+ C\ /c 2 is the (roughly constant) 
relativistic Alfven velocity, which allows us to use dt ~ — Va8 z . With these approximations we can reduce the system 
©-(113 to 



(d 2 -c 2 d 2 )B x 



E 



mn 



V A d z 



Bo + B x 



, 



(v z - V A )d z ( 7 v y ) - ±[v z (B + B x ) - V A B X ] = , 

TO 

{v z - V A )d z { 1 v z ) + ^-v y {B + B x ) = . 



m 



(17) 

(18) 
(19) 
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Since lu c 1 dt <C 1 and "fv y /c -C 1, Eq. (p~8|) can be used to obtain an approximate expression of v z , 



V A - 



Bo + B x 

which in turn can be used together with (|17p to express the wave equation as 

mn 2a B d z B x 



(d 2 -c 2 d 2 z )B x + ^^Vp. 



" [B 2 + 2B B x + (l-V 2 /c 2 )B 2 x f 2 



= 



By using d 2 — c 2 d 2 



where 



(c + VA){dt + cd z )d z we obtain 

d t B x + V{B x )d z B x = , 



V(B X ) = I c - (c - V A ) 



B x 



Vi\Bl 
Bl 



-3/2' 



(20) 



(21) 



(22) 



(23) 



We note here that the velocity V(B X ) agrees with the velocity U{B X ) in Eq. ([15]) (used previously by Ref. [T(|) to 
first order in c/Ca- Expression l|22p is slightly more accurate, however, as c/Ca *C 1 has not been used in the later 
derivation. 

Making the transformation t — r , z = ( + J Q T V(B((, r'))rfr' we have dt — ■> d T — (V(B x )/R)0q and d z —* (l/i?)9^, 
where R = 1 + J Q 7 (dV(B x ) / 'dQdr' '. As a consequence the general solution of Eq. (f2"2"| can be written t) = B x ((), 

which formally shows that wave steepening continues until wave breaking occur, such that R — > and d z — > cxd at 
some point for a finite time. Thus, in order to include dispersive effects that will prevent wave breaking we will need a 
more accurate expression than Eq. (|20[) for v z . Dispersive effects are then included by keeping terms of a higher order 
in an w^c^-expansion. How to implement this is to some extent dependent on whether we study an electron-ion 
plasma (in which case it is the higher order expression for the ion polarization drift that first leads to wave dispersion), 
or whether we consider an electron-positron plasma (in which case electron and positron motion contributes to wave 
dispersion simultaneously). For the astrophysical applications to be considered below, an electron-ion plasma is more 
appropriate, and thus we will focus on this case in the remainder of the article, using index i to denote ion quantities 
below. The only consequence for the calculations that has been made above is that we can simplify the expression 
for the Alfven velocity slightly, such that Ca = c (^2s w pl u> c)~ 1 ^ 2 = (^o/Mo^i^o) 1 ^ 2 - Next, using Eq. (|2TJ)) as an 
approximation of v Z i in Eq. (fl9|) gives 



nit 



Bin 

Bo + B x 



d z B x 



Bq + 2B B X 



1 - 



B 2 



3/2 



(24) 



which can be inserted into Eq. (|18p . resulting in a more accurate expression for v z 



V 3 R 3 
v A D o 



V A B X 



B Q + B x (B + B x 



d 2 z B x 



B 



J3 , 



{o z B x y 



Bl + 2BqB x +(l-Ypj Bl\ \B 2 + 2B B X + (l - Ypj B% 



This corrected expression of v Z i combined with Eq. (I17|) results in a wave equation of the form 



d t B x + V{B x )d z B x 



Va (Va 



Bo 



B x \u>, 



Bl 



(Bo + B x ) 



d 2 z B x 
9 5 (B X ) 



Bo 



v 2 

l - ^4 



B, 



(d z B x 



9 7 (B X ) 



= 



where V(B X ) is given by (|23[) and the auxiliary function g(B x ) is defined by 



g{B x ) 



Bq + 2B B X 



Yl 

c 2 



B\ 



1/2 



(25) 



(26) 



(27) 



As we can see, Eq. (|26|) now describes the fully nonlinear evolution combined with dispersive effects. In the weakly 
nonlinear limit this reduces to the celebrated Korteweg de Vries (KdV) equation, 



B x 



V 2 

d t B x + V A d z B x + 3(c - V A )^d z B x + (c- V A )^d 3 z B x 

Do LO^ 







(28) 
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As is wellknown [25(, for vanishing boundary conditions the general solution to the KdV equation involves a train of 
solitons, where the steepening effects due to nonlinearity is balanced by the dispersive effects due to the last term 
of Eq. (|28[). Our main interest here is the steepening of an initially sinusoidal wave profile. A study of the KdV 
equation for such initial conditions shows that the wave steepening induced by the nonlinearity leads to a harmonic 
content. However, as we are interested in generation of large harmonic content, with wave steepening occurring on a 
comparatively fast time scale, we need to study the fully nonlinear dispersive Eq. (j26|) rather than Eq. (|28|) . 

An initially sinusoidal wave described by Eq. (|26p can at first be approximated with Eq. (|22|) if the initial wave 
frequency lo fulfills lj <C uj C i- As a consequence the sinusoidal profile will undergo wave steepening until dispersive 
effects becomes important leading to a saturation of the steepening and the approach to a steady state profile, as can 
be described by Eq. (j2l))) . Thus we are looking for solutions to Eq. (f2T))) that are periodic and static in a frame moving 
with the wave. The general methods for finding solutions of this type is described in some detail in Ref. [2(|. First 
we let dt — > —VAndz, where Vau is a constant representing the nonlinearly modified Alfven velocity. Integrating the 
resulting equation with respect to z, we get an equation that are analogous to a particle moving in a potential [2q . 
For appropriate choices of integration constants, the motion corresponds to a particle moving in a potential well, in 
which case we obtain solutions of the desired type, i.e. B x = B x (z — VAnt), where B x is a periodic function of the 
argument. The resulting differential equation for B x has to be solved numerically. The degree of wave steepening 
in the resulting wave profile depends on the amplitude of the oscillations as well as on u>/uj C i, parameters which are 
related to the integration constants. For given integration constants, the corresponding physical parameters can be 
calculated straightforwardly, and thereby solutions with parameters of physical interest can be found by some trial 
and error. Solutions of astrophysical relevance showing significant wave steepening corresponding to large harmonic 
content of the wave profile will be presented in the preceding section. 



III. ASTROPHYSICAL EXAMPLE 



In this section we are going to investigate the nonlinearly modified Alfven waves, and to what extent these can 
be observed. Based on the fact that electromagnetic signals are easier to detect than gravitational ones of the same 
power, we will be looking for a spectral signature of GW-generated Alfven waves. As is wellknown, the gravitational 
sources have high powers only up to frequencies of the order ~ 1kHz, which is much below the lowest observable 
radio frequencies [U, [22]. Thus for GW induced radio wave detection to be possible, we will need a mechanism to 
convert some fraction of the electromagnetic energy to higher frequencies. As we will demonstrate below, the processes 
studied in Section II can provide the basis for such mechanisms, and thus lead to a large increase of the frequency for 
a significant fraction of the electromagnetic spectrum. 

Below we will study a concrete example. As a source of gravitational radiation we consider a binary system. At 
least one of the objects should have a strong magnetic field (in order to make the Alfven phase velocity close to c), and 
the objects should be compact (as to make the gravitational wave frequency and amplitude before merging reasonably 
large). For definiteness we study a system consisting of two neutron stars of equal mass Mq, separated by a distance 
of 20 Rs, where R$ = 2GM Q /c 2 ss 3 km. Furthermore, the surface magnetic field of each neutron star is assumed to 
be 4 x 10 6 T. The surroundings of the binary system can loosely be divided into three regions (Fig. [TJ, depending on 
which physical mechanisms that is dominating. 

I II III 




FIG. 1: The neighborhood of the binary system is divided into three regions: region I (10-Rs — 307?s), the energy conversion 
zone; region II (307is — 3500-Rs), the nonlinearity enhancement zone and region III (3500i?s — 7000Rs), the wave steepening 
zone. 
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A. Energy conversion zone 



The interval 10 Rs — 30 Rs from the center of mass (CM) roughly constitutes region I, which is the region where 
most of the gravitational energy is gained by the Alfven-wave. Using a Newtonian approximation, with d = aRg and 
r = (3Rs, h is straightforward to show that \h\ ~ (2a0)~ 1 , where d is the separation distance between the binary 
objects and r is the observation distance from the center of mass of the system. In order to obtain an estimate of 
the amplitude of the generated EM-wave, we note that in the near zone (i.e. where the magnetic pulsar magnetic 
field decays as a dipole) the plasma density is low, and we expect that for our purposes here, we can approximate the 
medium as vacuum. A calculation of the electromagnetic amplitude generated by the GW during such circumstances 
has been done by Ref. [|| . Combining the above expression for the gravitational wave amplitude at given distances 
with Eq. (14) of Ref. 0, the GW-induced magnetic field amplitude SB can thus be estimated. The result is 

^ ~ 1.8 x 10~ 4 , (29) 

DO 

at the of end region I. 



B. Nonlinearity enhancement zone 

In region II (approximate interval 30 Rs — 3500 Rs from the CM) the ratio determining the degree of nonlinearity, 
SB /Bo, is increasing quadratically . The reason is that SB suffers spherical attenuation (due to the high frequency) 
whereas Bo decays as a dipole field until the light cylinder is reached. The outer bound of region II, with a light 
cylinder at 3500i?s, is somewhat artificially chosen, and correspond to a pulsar period of 35 ms. Note, however, that 
the nonlinearity enhancement mechanism is sufficient to reach the strongly nonlinear regime (i.e. SB /Bo ~ 1) within 
region II even for faster pulsars. 



C. Wave steepening zone 

In region III (approximate interval 3500 Rs — 7000 Rs), the amplitude of the Alfven waves is strongly nonlinear, 
and thus pronounced wave steepening will occur here, as described by Eq. (|28|) . The characteristic propagation length 
L p for large steepening to occur is of the order 

where SVa is the velocity difference between the maximum and minimum wave velocity within the amplitude wave 
profile with dispersive effects neglected. At the beginning of region III we have Bo — 3 x 10 _3 T. Assuming a 
moderate plasma density, of the order no — 10 10 m~ 3 , gives an Alfven velocity Va — 0.9c. Furthermore, a relatively 
strong nonlinearity with SB /Bo = 0.25, consistent with the estimate made in section III B, correspond to SVa — 0.2c, 
which combined with a wavelength A ~ 1.7 x 10 6 m result in a steepening distance located reasonably well within 
region III. 

Due to the nonlinear velocity dependence of the amplitude, wave steepening continues until the dispersive effects 
become important. As the initial GW-generated angular frequency is uj « 1.1 x 10 3 s _1 , whereas Lu ci w 2.7 x 10 5 s _1 
in our case, and the amplitude is strongly nonlinear, the steady state profiles are dramatically changed as compared 
to the initially sinusoidal wave profiles. The steady state profile derived numerically from Eq. (f28|) corresponding 
to ui/LU C i — 1/240 and SB/B ~ 0.25 is shown in Fig. O Since, the temporal and spatial derivatives are orders of 
magnitudes larger than the corresponding values for the initial profile, the spectral content has changed dramatically 
from the initially quasi-monochromatic wave. The spectrum corresponding to the same data as in Fig. is shown 
in Fig. [3] There it is seen that the highest harmonic generation of Alfven waves may approach the 100 kHz range. 
Unfortunately, this is still not enough to reach the radio window for earth based antennas, where the lower cut-off lies 
around 20MHz. However, proposals such as the Astronomical Low Frequency Array (ALFA) [2l[ and more recently 
the Solar Imaging Radio Array (SIRA) [22| show that satellite based antennas can provide an opening for observation, 
since such arrays are planned for a sensitivity down to frequencies of the order 30kHz. As shown in Fig. [3l the spectral 
content exceeding this frequency is of the order 10% of the total wave energy density. Thus the astrophysical processes 
outlined here in section III opens for the possibility to correlate and interpret GW-observations made by LIGO or 
LISA with radio wave observations made with SIRA. 
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FIG. 2: The wave profile for a moderately nonlinear amplitude (SB/Bq « 0.25) as a function of the normalized distance 
z n = z luci/Va. The wavelength X n ~ 1500 correspond to a ratio cu/cud « 1/240. 
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FIG. 3: The normalized energy density (total energy density = 1) content W n of the wave profile as a function of the harmonic 
number n corresponding to the wave profile in Fig 2. The vertical line around n w f70 separates the harmonics exceeding the 
lower sensitivity bound of SIRA from the frequencies that are too low to be detected. For our given example, the amount of 
energy that can be detected is of the order of 10% of the total wave energy. 

IV. SUMMARY AND DISCUSSION 

We have considered the coupling of GW:s and Alfven waves in a magnetized plasma. The dependence of the 
Alfven velocity on the magnetic field amplitude makes the evolution equation for the Alfven waves strongly nonlinear. 
As described in section II, GW-generated Alfven waves are subject to nonlinear wave steepening, which is directly 
associated with harmonic generation. In order to find a saturation mechanism for the steepening process, we have 
included dispersive effects in our model by solving the momentum equation to third order in an <9 t /o; c -expansion. As 
demonstrated in section IIIc, the high amplitude steady state solutions for the Alfven waves show strong harmonic 
generation, which for sufficient amplitude does not saturate until the frequencies of the harmonics approach the 
cyclotron frequency of the ions. 

As a specific example we have considered a case where the GW:s associated with a neutron star-pulsar merging 
generate strongly nonlinear Alfven waves, resulting in high harmonic generation of the Alfven waves. As shown, the 
original GW-frequencies of the order ~ 200Hz may generate signals with harmonic number n ~ 200 or even larger. 
Unfortunately this is not sufficient to reach the radio window for earth based detection, but proposals for satellite based 
radio observations [22{ . that can detect radio signals down to ~ 30kHz could make observations possible. This opens for 
the future possibility to correlate GW-observations made by for example LISA [l[ with radio wave observations made 
with SIRA [22j |. to make detailed comparison with theories for the coupled electromagnetic-gravitational evolution. 
The process described here can be even more effective than in our example, if pulsars with stronger magnetic fields 
are considered. 

Finally we should point out that for the Alfven waves generated in our example to travel through interstellar 
distances, they must be more or less decoupled from the medium to become ordinary radio waves. For this to happen, 
the density should fall off with distance sufficiently fast after the wave steepening region in the previous example, such 
that the plasma frequency falls beyond the highest parts of the wave frequency before the cyclotron frequency does 
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so. In our example this is fulfilled, since the highest harmonic frequencies can exceed the plasma frequency already 
within the wave steepening region. However, for an accretion disc that is still thick at a distance when the magnetic 
field falls below the wave frequency, cyclotron damping could limit the possibilities for Radio wave detection to a 
significant extent. 
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